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•  Finite  element  models  are  developed  to  evaluate  the  effectiveness  of  the  pinch-torsion  test  in  creating  ISCrs  in  Li-ion  cell. 

•  The  torsion  component  can  trigger  ISCr  at  a  lower  load  with  smaller  short  spot  size. 

•  This  method  can  distinguish  commercial  cell  safety  performance  in  their  entire  range  of  SOC. 
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Recently  a  pinch-torsion  test  is  developed  for  safety  testing  of  Li-ion  batteries.  It  has  been  demon¬ 
strated  that  this  test  can  generate  small  internal  short-circuit  spots  in  the  separator  in  a  controllable 
and  repeatable  manner.  In  the  current  research,  the  failure  mechanism  is  examined  by  numerical 
simulations  and  comparisons  to  experimental  observations.  Finite  element  models  are  developed  to 
evaluate  the  deformation  of  the  separators  under  both  pure  pinch  and  pinch— torsion  loading  con¬ 
ditions.  It  is  discovered  that  the  addition  of  the  torsion  component  significantly  increased  the 
maximum  first  principal  strain,  which  is  believed  to  induce  the  internal  short  circuit.  In  addition,  the 
applied  load  in  the  pinch-torsion  test  is  significantly  less  than  in  the  pure  pinch  test,  thus  dramat¬ 
ically  improving  the  applicability  of  this  method  to  ultra-thick  batteries  which  otherwise  require 
heavy  load  in  excess  of  machine  capability.  It  is  further  found  that  the  separator  failure  is  achieved  in 
the  early  stage  of  torsion  (within  a  few  degree  of  rotation).  Effect  of  coefficient  of  friction  on  the 
maximum  first  principal  strain  is  also  examined. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  batteries  are  becoming  a  primary  power  source  in 
our  daily  lives  through  electronic  devices  (cell  phones,  tablets  and 
laptops)  and  transportations  (hybrid  and  electric  vehicles,  and 
airplanes).  These  applications  are  demanding  more  power  output, 
higher  power  density  and  lower  cost,  sometimes  at  the  expense  of 
safety.  Safety  hazard  related  issues  of  Li-ion  cells  have  been  well 
documented  [1,2].  Many  major  field  incidents  are  caused  by 
externally  or  internally  (e.g.,  manufacturing  defects,  mechanical 
abuse,  usage  abuse,  etc.)  induced  short  circuits,  which  can 
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potentially  release  the  high  energy  stored  in  the  battery  in  very 
short  time  locally  and  trigger  chemical  chain  reactions  releasing  a 
massive  amount  of  heat.  If  the  battery  is  not  well  designed  such 
that  the  heat  cannot  be  conducted  away  quickly,  thermal  runaway 
could  happen  and  lead  to  fires  and  explosions  in  some  extreme 
cases  [3,4].  Among  these  safety  concerns,  internal  short  circuit 
(ISCr)  under  no  obvious  abuse  or  external  triggers  are  less  under¬ 
stood  and  very  difficult  to  reproduce  experimentally.  Therefore, 
evaluating  the  thermal  runaway  risk  of  Li-ion  batteries  by  experi¬ 
mentally  creating  ISCr  in  a  controllable  and  predictable  manner  has 
brought  broad  interests  to  the  field. 

In  order  to  evaluate  the  risk  of  thermal  runaway,  many  tests 
have  been  conducted  to  simulate  ISCr  event  via  internal  defect 
initiation,  including  forced  ISCr  test  by  the  Battery  Association  of 
Japan  [5],  instrumented  indentation  [6,7]  and  nail  penetration  [8]. 
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Nomenclature 

a 

contact  radius 

m 

distance  between  the  maximum  tangential  traction 
location  and  the  contact  center 

Fo 

the  applied  load  on  the  indenter  in  pinch  only  test 

h 

the  thickness  of  cathode,  separator,  and  anode 
layers  in  the  unit  cell 

ISCr 

internal  short  circuit 

Mz 

the  twist  moment  applied  on  the  indenter 

P(r) 

normal  contact  pressure  distribution  in  the  contact 
area 

<7  (r) 

tangential  traction  in  the  contact  area 

r 

distance  to  the  contact  center  in  the  contact  area 

u i,  u2y  u3 

displacements  respectively  in  x,  y,  z  direction  in 

Fig.  2 

u4,  u5,  u6 

i  rotations  respectively  along  x,  y,  z  direction  in  Fig.  2 

£ 0 

the  maximum  first  principle  strain  in  pinch  only  test 
with  the  load  F0 

£\ 

the  first  principal  strain 

^Imax 

the  maximum  first  principal  strain 

£\\\ 

the  third  principal  strain 

£Y 

yield  strain 

e 

twist  angle 

friction  coefficient 

<7Y 

yield  strength 

On 

ultimate  strength 

We  recently  reported  an  improved  pinch  test  method  [9,10]  that 
could  reproducibly  create  controllable  ISCr  in  a  cell  separator  where 
the  size  of  the  ISCr  spots  depends  on  the  loading  speed,  pinch  ball 
or  indenter  diameter,  and  stroke  return-voltage.  A  further  devel¬ 
opment  of  this  pinch  test  method  added  a  torsional  loading 
component,  which  exhibited  improvement  in  the  effectiveness  of 
creating  the  ISCr  [11].  It  was  demonstrated  in  two  different  com¬ 
mercial  Li-ion  cells  that  the  torsion  facilitated  the  occurrence  of 
ISCr  with  lower  axial  load  and  smaller  ISCr  spot  size.  This  method  is 
thus  potentially  applicable  to  very  thick  batteries,  for  which  the 
critical  loads  under  pure  pinch  tests  are  too  large  for  the  standard 
loading  apparatus. 

In  order  to  quantitatively  relate  the  pinch-torsion  tests  to  the 
thermal  runaway  failure  in  batteries,  it  is  of  critical  importance  to 
understand  the  deformation  and  failure  mechanisms  under  such 
loading  conditions.  For  instance,  what  is  the  dependence  of  the 
critical  normal/twist  load  on  the  indenter  radius,  battery  thickness, 
number  of  electrode/separator  layers,  and  indenter/battery  friction 
condition,  among  many  other  factors?  Imagine  an  application  of 
successive  pinch  and  torsion  loads,  and  one  can  design  the  loading 
pattern/history  to  conveniently  identify  the  failure  initiation  and 
optimize  the  load  magnitude  to  allow  a  portable  evaluation. 
Motivated  by  the  above  considerations,  the  present  paper  attempts 
to  develop  a  computational  model  for  both  test  methods  (pinch- 
only  and  pinch-torsion  test  systems)  and  make  a  systematic 
investigation  on  the  deformation  mechanisms  to  provide  insight 
for  the  battery  safety  assessment.  In  this  paper,  finite  element 
method  (FEM)  is  used  to  simulate  the  strain  field  and  to  predict  the 
defect  initiation,  which  was  then  compared  with  experiment  re¬ 
sults.  The  effect  of  surface  condition  was  also  discussed  to  optimize 
the  test  design.  Moreover,  the  deformation  mechanisms  were 
rationalized  by  an  analytical  stick-slip  model  when  both  contact 
and  torsion  were  applied. 


Fig.  1.  Optical  graph  of  an  anode  layer  after  pinch  testing  showing  most  of  the  graphite 
coating  in  the  indented  region  was  removed  during  the  test  (the  gray  zone  is  graphite 
coating). 


2.  Numerical  model 

The  repeating  functional  unit  of  a  dry  Li-ion  cell  contains  an 
anode  layer,  a  cathode  layer,  and  a  separator  layer.  Since  electro¬ 
lytes  used  in  most  Li-ion  cells  are  in  liquid  phase,  they  are  not 
included  in  the  finite  element  models.  In  addition,  the  active 
coatings  (graphite  and  lithium  salts)  are  generally  loosely  bonded 
powders  on  the  electrodes.  During  slow  pinch-torsion  tests,  the 
coating  materials  were  delaminated,  worn  into  small  pieces  and 
displaced  from  the  highly  stressed  zone  near  the  indenter  tip.  In 
fact,  our  experimental  observation  revealed  that  the  active  coating 
materials  were  pushed  away  from  the  contact  area  after  testing 
(Fig.  1).  Therefore,  the  coatings  structurally  bear  little  load  espe¬ 
cially  comparing  with  the  other  components  and  will  be  simply 
excluded  in  our  model.  The  stress  equilibrium  equation  is  da y/ 
6x;  +  bj  =  0  in  which  bj  is  the  body  force  and  is  ignored  in  our  model. 
The  strain  ey  are  the  first  derivative  of  the  displacement  i.e. 
eij  =  ViidUildXj  +  dUj/dXi )  [13]  and  stress-strain  relationship  is  given 
by  the  constitutive  laws  which  will  be  discussed  later.  The 
boundary  value  problem  is  then  resolved  with  commercial  finite 
element  package  ABAQUS  (3DS  SIMULIA). 

Fig.  2  depicts  a  three-layer  (anode-separator-cathode)  unit  cell 
model  in  the  FEM  simulation  in  this  study.  The  material  of  the  top 
anode  layer  is  Cu,  the  separator  layer  in  the  middle  is  high-density 
polyethylene  (FIDPE)  and  the  bottom  cathode  layer  is  AL  The 
thickness  of  each  layer  h  is  0.02  mm.  The  two  indenters  in  the  pinch 
and  torsion  test  are  modeled  as  two  rigid  spheres  and  the  indenter 
radius  is  12.7  mm.  The  lateral  dimension  of  the  unit  cell  ( - 10  mm) 
is  chosen  to  be  ten  times  larger  than  the  estimated  contact  radius 


Fig.  2.  Schematics  of  the  cross  section  view  of  a  three-layer  battery  unit  cell  system 
under  the  pinch  tests  between  two  spherical  indenters. 
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Table  1 

Material  properties  of  Al  and  Cu  foils  measured  from  uniaxial  tension  tests  [12]. 


Aluminum  foil 

Copper  foil 

Yield  stress,  aY  (MPa) 

180 

210 

Ultimate  stress,  <ru  (MPa) 

195 

230 

Yield  strain,  ey 

0.017 

0.022 

Elongation  (%) 

2.7 

3.5 

size  (-1  mm)  to  eliminate  the  boundary  effects.  The  bottom 
indenter  is  fixed  (i.e.,  u2-  =  0,  i  =  1,  2... 6).  The  top  indenter  is  first 
applied  with  a  translation  123  =  0.9 h  vertical  to  the  unit  cell  top 
surface  and  then  rotated  with  displacement  boundary  condition 
Us  =  6  to  represent  the  twisting  in  the  experiment.  All  other  degree 
of  freedoms  of  the  top  indenter  are  zero. 

Both  Al  and  Cu  materials  are  modeled  as  elastic-plastic  solids 
with  isotropic  hardening  laws.  Materials  parameters  used  in  the 
finite  element  model  are  listed  in  Table  1.  The  HDPE  material  is 
modeled  by  the  hyperelastic  Neo-Hookean  material  model,  where 
the  “strain  energy  potential”  U  is  based  on  strain  invariants.  The 
form  of  the  Neo-Hookean  strain  energy  potential  is 

U  =  c,o(h-3)+±(j«-l)2  <» 

where  7^  is  the  first  deviatoric  strain  invariant  and  Jel  is  the  elastic 
volume  ratio,  C10  and  D\  are  the  model  parameters  determined 
from  the  HDPE  material  stress— strain  curve  measured  from  uni¬ 
axial  tensile  test  [13].  HDPE  materials  show  various  behavior  in 
uniaxial  tensile  test  e.g.  monotonic  hardening  [14]  and  softening 
effect  where  the  stress  decreases  with  increasing  strain  after  yield 
due  to  the  reorientation  phenomenon  [15].  Representative  stress- 
strain  curves  of  HDPE  are  shown  in  Fig.  3.  In  our  analysis,  discrete 
data  points  from  uniaxial  tension  test  are  input  in  the  material 
definition  section  in  ABAQUS. 

In  order  to  compare  with  the  experimental  procedure  and 
research  on  the  twisting  effect  to  the  defect  initiation  in  the  HDPE 
layer,  the  pinch  and  torsion  tests  are  modeled  in  two  successive 
steps  in  the  FEM  simulation.  In  the  first  step,  the  top  indenter 
vertically  pinches  on  the  top  surface  of  the  unit  cell.  This  is  an 
axisymmetric  problem  and  the  analysis  cost  is  substantially 
reduced  by  using  axisymmetric  elements  compared  with  a  fully 
three-dimensional  model.  In  this  simulation,  the  displacement  of 


Fig.  3.  Representative  stress-strain  curve  of  the  HDPE  material. 


the  top  indenter  is  0.018  mm  which  is  0.9h.  Large  deformation  (i.e., 
geometric  nonlinearity)  is  used  in  the  simulation.  In  the  second 
twisting  step,  torsion  along  z  axis  is  not  symmetric  with  z  axis  (axis 
of  symmetry  in  pinch  test)  and  thus  tzq  is  nonzero  in  the  battery. 
Therefore  a  full  three-dimensional  model  is  used  to  simulate  the 
indenter  twist  in  the  experiment.  The  displacement  and  stress 
fields  result  from  step  one  are  first  imported  into  the  three- 
dimensional  model,  and  the  top  indenter  is  then  applied  with  a 
maximum  rotation  angle  of  12°.  The  rigid  body  motion  of  the 
specimen  is  removed  by  fixing  the  transverse  translation  of  the  unit 
cell  center  line,  i.e.  u2  =  0,2  =  1,  2,  3.  The  Coulomb  friction  model  is 
used  between  the  indenter  and  the  specimen,  and  the  effect  of 
friction  coefficient  will  be  studied  in  details. 

Since  this  paper  focus  on  the  initiation  of  ISCr  created  by  the 
pinch  and  torsion  tests,  the  surface  interactions  within  the  unit  cell 
are  not  considered.  Therefore,  we  tie  each  contact  surfaces  within 
the  unit  cell,  i.e.  surfaces  between  the  Al  and  separator  layers  and 
surfaces  between  Cu  and  separator  layers  to  reduce  the  computa¬ 
tion  cost  from  the  contact  iterations  and  numerical  divergence  in 
ABAQUS. 

3.  FEM  results  and  discussion 

3.1.  Pinch  test 

In  the  pinch  and  torsion  experiments  [11],  ISCr  in  the  Li-ion 
battery  is  intentionally  created  by  breaking  the  separator  layer  to 
trigger  the  chemical  reactions.  The  cell  was  charged  to  a  set  level 
and  the  open  voltage  between  the  anode  and  the  cathode  is  closely 
monitored  to  ensure  that  the  indenter  will  be  retreated  immedi¬ 
ately  whenever  a  selected  voltage  drop  (e.g.  0.1  V)  is  measured. 
There  was  no  additional  charge/discharge  occurred  during  the 
mechanical  test.  There  is  some  volume  change  during  charge/ 
discharge,  which  adds  mechanical  stress.  However,  since  the  cur¬ 
rent  research  focus  on  externally  applied  mechanical  force  to 
induce  failure,  the  electrochemistry  induced  mechanical  behavior 
is  not  studied  in  this  paper.  The  failure  of  the  separator  is  due  to  the 
high  tensile  strain  caused  by  the  contact.  To  predict  the  ISCr  loca¬ 
tion,  the  tensile  zone  in  the  separator  layer  in  the  unit  cell  beneath 
the  spherical  contact  is  characterized  by  the  absolute  value  of  the 
first  principal  strain  over  the  third  principal  strain,  i.e.  \e\/em\,  and  is 
plotted  as  a  top  view  in  Fig.  4a.  Experimental  research  shows  that 
the  ductile  failure  process  in  HDPE  occurs  due  to  the  localized 
necking  after  excessive  elongation  at  an  intermediate  strain  rate,  a 
typical  rate  for  our  pinch-torsion  test  16].  When  \ei/sm\  >  1,  the 
location  is  under  tension  state  where  the  necking  is  mostly  likely  to 
initiate  (denoted  as  red  zones  in  Fig.  4a).  The  blue  zone  in  Fig.  4a 
undergoes  compressive  deformation.  It  is  shown  that  the  spherical 
indentation  results  into  a  local  tension  zone  in  the  separator  near 
the  indenter  but  outside  the  contact  area.  The  inner  radius  r  of  the 
tension  zone  is  about  1.2a  and  the  outer  radius  is  approximate  1.8a 
where  a  is  the  contact  radius.  The  distribution  of  the  first  principal 
strain  e\  in  the  separator  is  shown  in  Fig.  4b.  The  highest  value  of 
first  principal  strain  e\  in  separator  layer  exists  at  the  point  1.35a 
away  from  the  center  and  is  under  tension.  This  simulation  result 
qualitatively  agrees  with  our  experimental  examination  such  that 
holes  on  the  separators  usually  appeared  at  locations  some  distance 
away  from  the  center  of  the  impression  as  shown  in  Fig.  5.  After 
pinch  (Fig.  5a)  and  pinch  torsion  test  (Fig.  5b),  the  separators  layer 
in  the  unit  cell  are  observed  with  a  hole  outside  the  contact  area. 

3.2.  Pinch-torsion  test 

It  has  been  observed  that  some  initial  axial  load  needed  to  create 
ISCr  in  the  unit  cell  drops  significantly  if  a  small  amount  of  torsion 
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abs(LEl/LE3) 
(Avg:  50%) 

+  1.236e+01 
+3.000e+00 
+  2.833e+00 
+  2.667e+00 
+  2.500e+00 
—  +2.333e+00 

■-  +2.167e+00 
+  2.000e+00 

I+1.833e+00 
+  1.667e  +  00 
+  1.500e+00 
+  1.333e+00 
-M-  +1.167e+00 
+  1.000e+00 
+8.824e-04 


(a) 


LE,  Max.  Principal 
(Avg:  75°/o) 

+1.879e-02 

■  +1.723e-02 

■  +  1.566e-02 
—  +1.410e-02 

+  1.253e-02 
+  1.096e-02 
+9.397e-03 
+  7.831e-03 
+6.265e-03 
+4.699e-03 
, — ,  +3.132e-03 
^  +1.566e-03 

■  +7.542e-08 


(b) 


Fig.  4.  (a)  Top-view  contours  of  absolute  value  of  ei/em  in  the  separator  layer  in  pinch-only  test.  The  non-blue  colors  denote  the  region  under  tension,  (b)  Cross  section  view  of  the 
distribution  of  maximum  first  principal  strain  e\  after  an  indentation  of  0.018  mm  in  unit  cell.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred 
to  the  web  version  of  this  article.) 


is  applied  on  the  indenter  in  our  recently  developed  pinch-torsion 
test  method  11].  To  investigate  the  effect  of  torsional  component 
on  initiating  defect  during  pinch  test,  a  full  three-dimensional 
torsion  model  is  simulated.  As  stated  earlier,  we  first  import  data 
from  the  pure  pinch  test  and  then  conduct  a  static  analysis  by 
applying  a  rotation  on  the  indenter.  Fig.  6  shows  the  comparison  of 
the  strain  fields  in  the  separator  in  top  view  before  (Fig.  6a)  and 
after  (Fig.  6b)  the  torsion  is  applied.  The  friction  coefficients  be¬ 
tween  the  indenter  and  the  unit  cell  top  and  bottom  surfaces  are 
chosen  to  be  0.3.  After  the  torsion,  the  size  of  the  tension  zone  is  not 
changed  noticeably.  Flowever,  the  maximum  first  principal  strain  e\ 
in  the  tension  zone  increases  significantly  (by  77.7%)  after  the 
indenter  is  twisted  by  12°.  It  demonstrates  that  a  slight  twisting  of 
the  indenter  can  assist  the  defect  generation  which  is  consistent 
with  our  experimental  observation  [11  .In  order  to  understand  the 
effect  of  friction  coefficient  between  the  contact  surfaces,  p,  on  the 
pinch  and  torsion  test,  the  maximum  first  principal  strain  e\  is 
plotted  as  a  function  of  6  in  Fig.  7  for  three  different  friction  coef¬ 
ficient  values.  When  a  small  twist  is  applied  to  the  indenter,  e\  in¬ 
creases  almost  linearly  with  the  twist  angle.  The  slope  of  e\  and  6 
curve  gradually  decrease  with  increasing  twist  angle  until  e\  rea¬ 
ches  a  plateau  at  a  twist  angle  about  5°. 

The  above  phenomenon  can  be  explained  by  an  elastic  contact 
analysis  [17-19].  For  an  elastic  half-space  under  spherical  contact, 
the  normal  contact  pressure  distribution  in  the  contact  area  is  given 
by  the  Flertz  theory: 

P(r)  =  Po/l  -  Q2  (2) 

where  a  is  the  contact  radius,  po  =  (6PE*2/n3R2y^3  is  the  maximum 
contact  pressure  and  E*  is  the  reduced  modulus  [20].  The 


(a) 


(b) 


Fig.  5.  Optical  graphs  showing  broken  separators  after  (a)  a  pinch  test  and  (b)  a  pinch- 
torsion  test.  The  pointed  cracks  in  the  separator  layer  are  initiated  outside  the  spherical 
contact  surface  for  both  tests. 


360 


Y.  Xia  et  al.  /  Journal  of  Power  Sources  265  (2014)  356-362 


(a) 


LE,  Max.  Principal 
(Avg:  75%) 
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(b) 


Fig.  6.  Simulated  distribution  of  e\  in  separator  layer  (a)  before  twisting  (b)  after  twisting  with  angle  of  12°.  The  black  dash  lines  denote  the  tension  zone. 


maximum  contact  pressure  resides  on  the  center  of  the  contact 
area  and  monotonically  decreases  along  the  radial  direction.  If  a 
torque  Mz  is  applied,  the  tangential  traction  q(r)  in  the  contact  area 
is  distributed  as, 


<?(»') 


3  Mzr 
4na3 


(a2 


(3) 


where,  mathematically,  the  torsional  shear  traction  goes  to  infinity 
at  the  edge  of  the  contact  area  (i.e.  r  -►  a).  Such  a  solution  is  similar 
to  the  Mode-Ill  crack  tip  field  at  the  contact  edge  where  the  classic 
inverse-square-root  singularity  resides,  and  obviously  it  holds  only 
when  the  interface  is  perfectly  bonded.  According  to  the  Coulomb 
friction  law,  when  \q(r)\  <  p|p(r)|,  there  is  no  tangential  slip  be¬ 
tween  the  two  contact  surfaces,  and  when  |q(r)|  >  p|p(r)|,  there  will 
be  a  relative  slip  between  the  indenter  and  the  substrate.  The  shear 
stress  singularity  at  the  contact  edge  clearly  leads  to  the  stick-slip 
behavior  as  shown  in  Fig.  8a.  From  Eq.  (3),  after  a  torque  Mz(0)  is 
applied,  where  6  is  the  twist  angle,  the  edge  of  the  contact  surface 
will  start  to  slip  immediately  due  to  the  large  tangential  traction. 
However,  the  rest  of  the  region  in  the  contact  area  remains  bonded, 
making  the  indenter  contact  surface  partially  slip  on  the  solid.  For 
any  point  in  the  contact  area  which  has  a  distance  r  to  the  contact 
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Fig.  7.  The  maximum  first  principal  strain  e\  as  a  function  of  twist  angle  6  with 
different  friction  coefficients  /l.  The  strain  e\  during  the  torsion  test  is  normalized  by 
the  maximum  first  principal  strain  in  HDPE  when  the  pinch  test  ends  which  has 
negligible  dependence  on  the  friction  coefficient. 
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(b) 

Fig.  8.  (a)  Shear  traction  q(r/a)  distribution  when  the  indenter  is  under  partial  slip,  (b) 
The  strain  eze  distribution  along  the  radius  at  the  contact  surface  when  6  =  2.25°  and 
7.2°. 

center,  the  critical  twist  angle  needed  to  initiate  the  slide  is  related 
to  indent  pressure  as, 

q(r)  =  np{r)  =  |^(a2  -  r2)2  (4) 

In  this  study,  we  assume  the  static  friction  coefficient  is  the 
same  as  the  kinetic  friction  coefficient.  Therefore,  once  the  two 
contact  surfaces  change  state  from  stick  to  slip,  the  tangential 
traction  will  remain  as  a  constant  at  the  value  of  maximum  static 
shear  force,  i.e.  /rp(r).  The  largest  q(r)  and  thus  the  largest  shear 
stress  tZ0  in  the  contact  surface  always  occurs  at  the  place  where  it 
is  just  start  to  slip.  We  denote  the  distance  between  the  q(r)  max¬ 
ima  and  the  contact  center  as  c{6).  A  representative  profile  of  q(r/a) 
when  the  indenter  partially  slips  on  the  contact  surface  is  plotted  in 
Fig.  8a.  Where  q(r)  reaches  its  peak  at  c{6).  For  location  r  >  c{6)  in 
the  contact  area,  it  is  in  dynamic  slipping  while  location  remains 
static  at  r  <  c{6).  The  shear  stress  tzq  will  add  a  large  shear  strain  yzo 
and  thus  increase  the  maximum  e\,  corresponding  the  ascending 
part  in  Fig.  4.  The  stick  circle  eventually  shrinks  towards  the  contact 
center  (i.e.  c{6)  -►  0)  and  the  free  sliding  of  the  indenter  is  started 
with  a  constant  moment  Mz.  Therefore,  no  further  deformation 
takes  place  in  the  battery  system  and  the  maximum  first  principal 
strain  e\  arrives  in  the  plateau  shown  in  Fig.  7. 


The  above  analytic  solution  is  derived  from  elastic  half-spaces. 
For  the  thin  layer  under  a  pair  of  spherical  contacts,  the  stick-slip 
behavior  applies  qualitatively,  but  its  quantitative  feature  can 
only  be  obtained  from  finite  element  analysis.  The  strain  ezq  is 
plotted  along  the  radius  at  the  battery  contact  surface  with  twist 
angles:  2.25°  and  7.2°  respectively  in  Fig.  8b.  For  6  =  2.25,  the 
maximum  shear  strain  appears  at  c{6)  =  0.3a.  For  6  =  7.2°,  the 
indenter  has  fully  slipped  so  that  max  ezq  appears  in  the  vicinity  of 
contact.  Based  on  the  analysis  above,  it  is  confirmed  that  torsion  of 
the  indenter  will  increase  the  intensity  of  tension  inside  the  sample 
thus  will  increase  the  possibility  of  breaking  the  separator  layer. 
Flowever,  such  effect  saturates  beyond  a  critical  twist  angle  and  e\ 
remains  relatively  constant,  which  depends  on  the  friction  coeffi¬ 
cient  /JL. 

To  explore  e\  plateau  value  in  the  separator  layer,  the  influence  of 
friction  coefficient  /jl  between  the  surfaces  is  presented  in  Fig.  9.  It 
shows  that  the  plateau  monotonically  increases  with  the  increase 
of  friction  coefficient  /jl  because  from  Eq.  (3),  larger  friction  coeffi¬ 
cient  will  lead  to  larger  tangential  stress  and  thus  larger  principal 
tensile  strain.  In  other  words,  when  /jl  is  large,  a  larger  torque  is 
needed  to  produce  a  specific  twist  angle  for  the  indenter,  and  thus 
the  separator  layer  will  have  more  severe  deformation  to  assist  the 
defect  initiation.  In  conclusion,  torsion  component  facilitates  the 
failure  of  the  sample  by  increasing  tensile  deformation  in  the 
polymeric  separator. 

Our  recent  study  compared  the  pinch  and  pinch-torsion  tests  in 
terms  of  the  axial  load  needed  to  induce  ISCr  [11  ].  For  a  three-layer 
dry  cell  that  had  a  configuration  similar  to  Fig.  2,  it  was  discovered 
the  axial  load  generating  ISCr  decreased  from  approximately  630  N 
under  pure  pinch  loading  to  about  400  N  under  a  combined  pinch/ 
torsion  loading  -  a  reduction  of  36.5%  [11].  The  reduction  is 
explained  in  Fig.  10  where  a  fix  twist  angle  is  applied  to  different 
indenter  force  F  and  the  maximum  principle  strains  in  different 
cases  are  plotted.  We  also  make  a  pinch-only  test,  as  a  reference 
case,  by  applying  an  indentation  force  F0  on  the  unit  cell  and  read 
the  maximum  first  principle  strain  eq  without  any  twist  to  the 
indenter.  Fig.  10  shows  that  after  a  slight  twist  angle  2.7°  is  applied, 
one  would  need  less  indentation  force  (a  reduction  of  33%)  to  create 
same  amount  of  maximum  first  principal  strain  (i.e.  e\Ieo  =  1 )  in  the 
unit  cell.  It  is  in  line  with  what  is  observed  in  our  recent  study  11  ]. 
On  the  other  hand,  it  is  noticed  that  at  the  same  applied  load,  the 


Fig.  9.  The  maximum  first  principal  strain  e\  after  the  indenter  is  twisted  with  an  angle 
6  =  6  as  a  function  of  friction  coefficients  /jl. 
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Fig.  10.  The  maximum  first  principal  strain  e\  as  a  function  of  indentation  force  at 
different  fixed  twist  angles.  F0  and  e0,  respectively,  are  the  indentation  force  and  the 
maximum  first  principal  strain  when  no  twisting  is  applied.  The  experimental  data 
point  was  taken  from  Ref  [11  ],  which  represented  results  from  multiple  measurements. 
Since  those  samples  were  nominally  identical  they  failed  under  similar  conditions  and 
appeared  to  be  a  “single”  data  point  in  Fig.  10,  which  has  normalized  axes. 


maximum  first  principle  strain  inside  HDPE  increases  significantly 
with  6  from  1°  to  2°.  The  coincidence  of  data  for  angles  larger  than 
2°  shows  that  the  indenter  is  in  dynamic  slipping  already  after 
twisted  2°  at  small  applied  load  level  (i.e.  F/F0  is  less  than  0.9). 
When  indenting  with  a  large  load,  e.g.  F/Fo  ->  1,  higher  twisting 
angle  is  needed  to  make  the  contacting  surfaces  fully  slip.  And  the 
maximum  first  principle  strain  increases  by  17%  with  the  twisting 
angle  6  from  2°  to  5°. 

4.  Conclusions 

In  this  paper,  finite  element  analysis  is  used  to  simulate  defor¬ 
mation  of  Li-ion  battery  layers  under  pinch  and  pinch— torsion 
conditions.  The  simulation  discovery  is  in  agreement  with  our 
previous  experimental  results.  Several  conclusions  can  be  drawn  as 
follows: 

1.  Spherical  indentation  delivers  tensile  zone  inside  the  sample 
instead  on  the  surface  comparing  with  other  experiments  such 
as  flat  end  punch  and  bending  test,  making  it  an  attractive 
method  to  simulate  ISCr  in  Li-ion  batteries. 


2.  A  slight  twist  of  the  indenter  introduces  the  shear  strain,  which 
increases  the  maximum  first  principal  strain  in  tension  zone  and 
thus  facilitates  the  failure  of  the  polymer.  During  this  stage,  the 
indenter  and  the  contact  surface  partially  slip  with  co-existence 
of  a  slip  region  and  a  stick  region. 

3.  Beyond  a  critical  angle,  the  indenter  is  free  to  spin  relative  to  the 
battery  and  the  maximum  first  principal  strain  e\  takes  its  limit 
value,  which  highly  depends  on  the  friction  coefficient  /i. 
Therefore,  controlling  surface  condition  on  the  indenter  or 
battery  surface  could  be  utilized  to  control  ISCr  events. 
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